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Abstract 



In this letter we clarify the role of heat flux in the hydrodynamic balance 

equations in 2D quantum wells, facilitating the formulation of an Onsager 

relation within the framework of this theory. We find that the Onsager relation 

is satisfied within the framework of the 2D hydrodynamic balance equation 

transport theory at sufficiently high density. The condition of high density 

is consonant with the requirement of strong electron-electron interactions for 

the validity of our balance equation formulation. 
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The Lei-Ting balance equation transport theory has achieved much success in hot- 
electron transport of homogeneous semiconductors. This theory was subsequently general- 
ized to deal with weakly nonuniform, inhomogeneous systems 0. The resulting hydrody- 
namic balance equations consist of continuity equation, momentum balance equation, and 
energy balance equation. A salient feature of this hydrodynamic approach is that it in- 
cludes a microscopic description of scattering in the form of a frictional force function due 
to electron-impurity and electron-phonon scattering, as well as an electron energy loss rate 
function due to electron-phonon interaction. These hydrodynamic balance equations have 
recently been applied to device simulations and they have been further developed to 
include phonon-drag effects and applied to discussion of thermoelectric power under both 
linear [0 and nonlinear transport conditions. 

Until recently, the important issue concerning the capability of this theory to lead to 
the correct form of Onsager relations and/or how to determine Onsager relations within 
the framework of this theory has not been addressed. There is even some misunderstanding 
that the energy flux predicted by this theory is zero. It is well known that the Onsager rela- 
tion is a manifestation of microscopic irreversibility for any statistical system near thermal 
equilibrium. Therefore any properly formulated statistical physics model should satisfy this 
relation. It is very easy to verify this relation in the framework of Kubo linear response 
theory. Moreover, if one can determine the distribution function from the Boltzmann equa- 
tion, it is also straightforward to verify the Onsager relation by calculating the pertinent 
moments of the distribution function. However, for the traditional hydrodynamic model, 
||T0|^T1 which is derived from Boltzmann equation, verification has been elusive. In fact, in 



a very recent article [^], it has been argued that the Onsager relation breaks down in this 
model. 



Recently we clarified the role of heat flux in this theory, and, by introducing the 
fourth balance equation, ze., energy flux balance equation, we were able to show how to 
generate Onsager relations within the framework of this theory. Moreover, we closely checked 
the Onsager relation predicted by this theory for bulk semiconductors and found, that for 



any temperature, when electron density is sufficiently high, the balance equation theory 
satisfies Onsager relations exactly. To our knowledge, this is the first set of hydrodynamic 
equations which has been shown to obey Onsager relation exactly. 

The purpose of the present letter is to clarify the role of heat fiux in hydrodynamic 
balance equations in quantum wells (and other 2D systems), facilitating the formulation of 
an Onsager relation within the framework of this theory, and to verify the validity of this 
theory in regard to the Onsager relation. 



Following the procedure set forth in Ref . p3[ , the hydrodynamic balance equations which 



describe a weekly inhomogeneous electron system under the infiuence of an electron E and 
a small lattice temperature gradient VT in two dimensional (2D) quantum wells can be 
written as 

^ + V-(nv)=0, (1) 

|(J) + V-((J)v) = -lv« + i, (2) 
ot mm 

du 11 

— — h V • (J//) = V ■ Vu H — mv'^V ■ (rav) H — mnv ■ Vw^ — w — v ■ f , (3) 

which describe continuity, momentum balance and energy balance respectively. In these 
equations, the carrier drift velocity v, the electron temperature Tg, the average relative 
electron energy m, the carrier density n and the chemical potential /z, together with the 
particle fiux (J) and energy fiux (J//) are all field quantities weakly dependent on the 
spatial coordinate, such that their spatial gradients are small and therefore are retained 
only to first order in Eqs. (|l|-^). The energy fiux {iu) appearing in Eq. (|) is just the 
energy fiux predicted by hydrodynamic balance equation theory. It relation to the other 
field quantities in 2D case can be derived following Ref. and be written as 

(Ji^(R)) = 2m(R)v(R) + ^mn(R)t;2(R)v(R) . (4) 

Substituting this relation, together with the relation of density fiux 

(J(R)) = n(R)v(R) , (5) 



into Eqs. (Q) and one arrives at the original hydrodynamic balance equations p4| 



In order to obtain the Onsager relation in the framework of balance equation theory, we 
need first derive the energy-fiux balance equation [^. This in the case of 2D electrons in 
quantum wells, is given as 

^{Jh) + V ■ {A) = (B) + —euE + enE ■ vv + ^env'^E + ^v'^i - wv . (6) 

The expression of (B) is composed of two parts. One is due to collisions with impurities 
((Bj)), and the other is due to interaction with phonons ((Bp/i)). They are given by 

k + q/2 

(Bi) = 27mi ^ |M(q)nF(q, 2;o)P(ek+q - ^k) — ^^^^^(ek+q - ^k + q ■ v) 



kq 

X 



/( 



gk - /^ ^ 



/( 



k ~\~ 

+ 27rni V |M(q)nF(q, ^o)r(q ■ v h k • v— )(5(£k+q - ^k + q ■ v) 

mm 



kq 



/( 



f^kH 



Tin 



|M(q)nF(q,2;o)P(£:k+q + £:k)— 5(£k+q - £^k + q ■ v) 



kq 

X 



m 



(7) 



and 
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The tensor A can be expressed as 



{A) = -[S'(R) + + Smvv + -mnv'^w 



(9) 



with 

(10) 

In these equations = k'^/2m, and f{x) = l/(e^ + 1) represent the energy dispersion of 2D 
electrons and Fermi distribution function separately, and I stands for the unit tensor. 



The Onsager relation is concerned with the linear response of the particle current 
(J) and the heat flux (Jq) near thermal equilibrium, which flow as a result of generalized 
forces Xj on the system: 

{J) = L'^Xi + L^^X2 , (11) 
(Jq) = L^^Xi + L'^X, , (12) 

with Xi = —j^V^fi + ecj)) and X2 = V(l/T). The Onsager relation states that 

L'^ = . (13) 

The heat flux (Jq) relates to the energy flux in Eq. (^) through 

(Jq) = (Jh)-/x(J). (14) 

The fluxes (J) and (Jh) have already been defined by Eqs. (^ and (§). Our task is to express 
them in terms of linear response in the form of Eqs. ([Tl| ) and ([T^) . To this end, we consider 



electron transport in a quantum well which is grown along 2r-direction), in the presence 
of a small lattice temperature gradient and a small electric field along the x direction: 
VT = (VxT,0,0) and E = {Ex,0,0) respectively. Therefore Te = T, and v = {vx,0,0) is 
also small. The system is near equilibrium. Following the same steps as in Ref. we 
can derive the first relation (|TT]) from the momentum balance equation (0) and the second 
relation ([T2|) from the energy- flux balance equation (H), by treating Eqs. (0) and (H) to first 
order in the small quantities E^, V xT and v^, and have 
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In deriving these equations, we have used the relations 
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(20) 



with ( = p/T and Fi,{y) = /o°° x'^[exp(x — y) + 1] ^dx. Here, 

P = = ^ (21) 

is the resistivity and independent of {{Jx})- The expression of it can be found in Ref. 0. 
Further, 

l/r = ^^"^ (22) 

is also independent of Vx (therefore {Jfj))- Comparing Eq. (|l^) with Eq. ([T7|), we find that 
the condition under which the Onsager relation holds is given by 

-2 



7^-^ = 1 
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(23) 



We have closely examined Eq. (p3|) for a GaAs-based (10 nm well width along the z- 
direction) quantum well structure. (We have also considered other quantum wells of different 
well width and found that the tendency of / approaching to "1" is not very sensitive on the 
width of quantum well.) We only take the lowest subband occupation into account. Both p 
and (Bx) are composed of contributions due to electron-impurity, electron-LO-phonon, and 
electron-acoustic-phonon scatterings (with the electron-acoustic-phonon scatterings due to 
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longitudinal mode acoustic phonons via deformation potential and piezoelectric interac- 
tions, and transverse mode via piezoelectric interaction). We have examined each scattering 
contribution in detail to check Eq. p3|) separately for each interaction. It is clear that if 



L = —^TTT^ = 1 is satisfied for each interaction, we have — ^Pfrr^ = 1. Moreover, this 
procedure is advantageous in that each Jj is independent of impurity concentration and 
parameters of the electron-phonon interaction matrixes. 

The expressions for / obtained from the balance equations are given by 

(24) 
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due to electron-impurity scattering; and 
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due to electron-phonon scattering, for phonons of mode A. Ie-ph{^) is further composed 
of contributions due to electron-LO-phonon scattering, le-io'-, due to electron-longitudinal 
acoustic phonons by deformation potential coupling, ledu and by piezoelectric interaction, 
lepu and due to electron-transverse acoustic phonons by piezoelectric interaction, If,pt- In 
these equations, ^qa is the phonon frequency of wave vector Q = (q, gz) = {qx,%,<lz) 
and mode A; u{q), the electron-impurity interaction potential; and M(q, A), the electron- 
phonon interaction matrix element. -F(q, Zq) and I{iqz) are form factors of electron-impurity 
and electron-phonon interaction respectively, with zq standing for the position of impurity. 
n{x) = (e^^ — 1)~^ stands for the Bose distribution. 112 (q, ^) denotes the imaginary part of 
electron density- density correction function. n| is defined by 

n^2(q,^) = 2vrX:^k5(ek+,-ek + a;)[/(^)-/(^^^±^)] . (26) 

For the LO phonon, VLq^lq = = 35.4 meV, and the Frohlich matrix element is 
|M(Q,L0)|2 = e^{K^ - K-^)no/{2eoQ^) oc (Since the constants in the matrix 

elements cancel in Eq. (p5|) , in the following we only specify their relation to Q.) The 
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matrix element due to longitudinal deformation potential coupling is |M(Q,(i/)p oc Q, 
that due to longitudinal piezoelectric interaction is |M(Q,pZ)p oc isixflyflzf' I Q'^ ■, and for 
the two branches of independent transverse piezoelectric interaction: Z]j=i,2 l^(Q)P^j)P oc 
(O'^Q'y + + ^z^x ~ {^Ix^y^zY / Q"^) / ■ For acoustic phonons ^qa can be written as VsQ, 
with the longitudinal sound speed Vs being 5.29x10^ m/s, and the transverse sound speed 
being 2.48x10^ m/s. The effective mass of electron is 0.07me, with nie denoting the free 
electron mass. 

We present the results of our numerical calculations in Fig. 1 to Fig. 4, where contribu- 
tions to / due to the various interactions discussed above are plotted against electron density 
for several different temperatures (the results of I^pi are similar to those of /gpt, so we only 
plot one as representative). As it is generally believed that the contribution of acoustic 
phonons is important only at low temperature, while the contribution of LO phonons is 
dominant at high temperature, our temperatures are chosen as 10, 20, and 40 K for the 
former, and 50, 300, 500, and 1000 K for the letter. Impurity scattering is present at all 
temperature, so we take T =10, 50, 100, 300, and 1000 K in Fig. 1. From these figures it 
is evident that, for any temperature, when electron density is sufficiently high I is exactly 
unity, indicating that the Onsager relation holds. This is consistent with the understanding 
that the Lei-Ting balance equation theory holds only for strong electron-electron interaction 
21,11. 
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FIGURES 

FIG. 1. / due to electron-impurity scattering is plotted as a function of electron density for 

several different temperatures 

FIG. 2. / due to electron-LO-phonon scattering is plotted as a function of electron density for 
several different temperatures 

FIG. 3. I due to electron-longitudinal acoustic-phonon scattering via deformation potential 
coupling is plotted as a function of electron density for several different temperatures 

FIG. 4. I due to electron-transverse acoustic-phonon scattering via piezoelectric interaction is 
plotted as a function of electron density for several different temperatures 



11 



